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Abstract. The effective dynamics of solitons for the generalized nonlinear 
Schrodinger equation in a random potential is rigorously studied. It is shown 
that when the external potential varies slowly in space compared to the size of 
the soliton, the dynamics of the center of the soliton is almost surely described 
by Hamilton's equations for a classical particle in the random potential, plus 
error terms due to radiation damping. Furthermore, a limit theorem for the 
dynamics of the center of mass of the soliton in the weak-coupling and space- 
adiabatic limit is proven in two and higher dimensions: Under certain mixing 
hypotheses for the potential, the momentum of the center of mass of the soliton 
converges in law to a diffusion process on a sphere of constant momentum. 
Moreover, in three and higher dimensions, the trajectory of the center of mass 
of the soliton converges to a spatial Brownian motion. 



1. Introduction 

1.1. Overview of earlier results and heuristic discussion. In the last few 
years, there has been substantial progress in rigorously understanding solitary 
wave dynamics for the nonlinear Schrodinger equation in slowly varying potentials 
(or in the presence of small rough perturbations), see PQ-[9]. The basic picture is 
that in the space- adiabatic limit, or in the presence of small perturbations, the 
long-time dynamics of the center of the soliton is described by Hamilton's (or 
Newton's) equations in an effective potential that corresponds to the restriction 
of the external potential to the soliton manifold, plus error terms due to radiation 
damping. 

In the above cited work, a soliton behaves like a classical point particle over 
certain scales. The main question that we address is whether this analogy between 
solitons and point particles over a certain temporal and spatial scale still holds 
when the soliton is moving in the field of a random potential. We show that 
the answer to the above question is affirmative. Furthermore, we have strong 
enough explicit control over the long-time dynamics of the soliton in the random 
potential that allows us to look at the limiting dynamics over large-distances and 
long-times. In a certain weak-coupling /space-adiabatic limit, and under some 
mixing hypotheses for the random potential, we show that the long-time, large- 
distance behavior of the soliton is described by momentum diffusion in N > 2 
dimensions: The soliton center of mass undergoes Brownian motion on the energy 
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sphere of constant momentum. This is analogous to the long-time/large-distance 
behavior of the a classical particle in a random potential, [2] - [16]. Moreover, 
in dimensions N > 3, the long-time limit of a momentum diffusion is a spatial 
Brownian motion, [T7j. We also show this spatial diffusive behavior for the soliton 
center of mass. We note that the weak-coupling and space-adiabatic limit are 
taken simultaneously. This is more difficult than taking the semi-classical limit 
first, and then the weak-coupling limit. 

This diffusive motion of the center of mass has been observed numerically in 
[TUt [TT] for the NLS equation with power nonlinearity and white multiplicative 
noise; see also [12] for a discussion of a concrete experiment where such a dynamics 
may be observed. For the derivation of the NLS equation in the mean-field limit 
of interacting bosons in a random potential, we refer the reader to |13j . 

We note that a somewhat related problem arises in the semi-classical limit of 
the dynamics of a quantum particle in a random field. It is shown by Erdos, 
Salmhofer and Yau in [18] that the semi-classical/weak-coupling limit of the dy- 
namics of a quantum particle in a random potential displays a diffusive behavior 
of the energy density of solutions of the linear Schrodinger equation. Another 
related problem is that of localization for the nonlinear Schrodinger equation, 
which has been studied recently by Bourgain and Wang in [T5] . 

1.2. Description of the problem. In what follows, we consider the probability 
triple (Q,J-,¥), such that the probability space Q has a generic point u and is 
endowed with measure //. For a measurable and integrable function / on Q, we 
define the expectation value of / as E(/) := j f{uj)n{duj). 

In this paper, we study the long-time dynamics of solitary wave solutions for 
the generalized nonlinear Schrodinger equation in a random potential 

(1) id t i/;(x, t) = (- A + \V h (x, t; w))^(x, t) - /(^(x, *)), 

where uo G Q, x G R denotes a point in the configuration space, t G M is time, 
dt = J^, A = YujLi afp' the iV-dimensional Laplacian, Vh is the random potential, 
which is a measurable real function on R N xix VI — > K satisfying 

V h (x,t;w) = V(hx,t;w), h G (0,1], 

A G [0, 1] is a coupling constant, and the nonlinearity / is a mapping on complex 
Sobolev spaces, 

/ : i? 1 (E 7V ; C) -> H-\R N ;C), 

with /(0) = 0, and f(ip) = f(4>), where 7 denotes complex conjugation. 
Typical focusing nonlinearities are local ones 



(2) 



When A = 0, the NLS equation with nonlinearity as given above admits 
solitary wave solutions, which are stable spherically symmetric and positive so- 
lutions 



(3) V*(x,t) := e i( 5 v -< x - a )+T) 



x — a 
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where a := (a, v, 7, fi), a = vt + a , 7 = fit + ^j-t + 70, with 70 G [0, 27r), 
a ,v G M , /i G R + , and 77^ is a positive solution of the nonlinear eigenvalue 
problem 

(4) (-A + n)^ - f( Vll ) = 0. 

The solution ([3]) stands for a soliton traveling with velocity v, center of mass 
position a, and phase ~v • (x — a) + 7. The size of the soliton is oc /i -1 ^ 2 , in the 
sense that 77^ ~ e^^"*" as ||x|| — > 00, see for example [27] and [28] . 

We consider in this paper an external potential whose realizations are P-a.s. 
varying slowly in space compared to the size of the soliton, i.e., 

ess — sup a;gCj igRj xgR iv —= <^ 1, 

which corresponds to the space-adiabatic limit (h <C 1) if we set the size of the 
soliton to 0(1). Furthermore, we assume that the realizations of the external 
potential satisfy 

(5) V h G W 1,00 (R, C 2 {R N ) n W 2,co (R N )) P-a.s., 
i.e., there exists HcO with /x(fi) = 1, such that 

V h , d t V h G C 2 (R N ) n H /2 '°°(R Ar )) for w gH. 

We are interested in the dynamics of a single soliton in the external potential. 
We assume that there exists <7o = (a , v , 70, /io) such that 

(6) \m=o)- VtTO \\ H i<c h, 

for some positive constant Co, where h appears in ([1]). Theorem [1] below is a 
rigorous result in the space-adiabatic limit for the long-time dynamics of the soli- 
ton for the special class of nonlinearities and potential discussed above. A more 
general result, Theorem HI is stated in Section [2] after listing general assumptions, 
and is proven in Section |H 

Theorem 1. Consider the NLS equation (TJJ) with nonlinearity f given by (TJ|) and 
potential Vh satisfying (TJP , and suppose the initial condition satisfies Then 
there exists h > and positive constants C and C, such that, for all h G (0, ho) 
and for any fixed e G (0, 1), we have that 

_sup E[\\^- Va{t) \\ H i]<Ch 1 -^ 

te[0,Ce|log/i|/A/i) 



uniformly in X 6 (h 1 e , 1], and the parameters a (t) = (a(i), v(t), 7 (t), satisfy 
the system of effective equations 

dtSL = v + 0(/i 2 " e ) 

d 4 v = -2AVV r fc (a, t; w) + 0(/i 2 ~ e ) 

d t7 = // + J||v|| 2 - 7 ft (a, t; u) + 0(/i 2 ^) 

9 t /i = 0(/i 2 - e ), 

with initial condition satisfying 

||a(0) - ao||, ||v(0) - v ||, | 7 (0) - 7 o|, 1^(0) - /i | = 0(/i). 

We now discuss a limit theorem for the dynamics of the center of mass of 
the soliton moving in a time-independent and strongly mixing random potential. 
Suppose that 

(7) V h (x,t;w) = V(hx.\u), u e SI. 
Performing the simple scaling 

a := ha 
v := v 

t := ht, 

it follows from Theorem [1] that, for t £ [0, Ce| log /i|/A), 

(8) c\a = v + 0(h 2 -') 

(9) dfv = -2AVF(a;cu) + 0(/i 1 ~ e ), 
with initial condition 

(10) ||a(0) - fcaoll = 0(h 2 ), ||v(0) - v || = 0(/i), 

which are, up to small remainder terms, the equations of motion of a classical 
particle in a random potential. This motivates introducing an auxiliary stochastic 
process (a(t), v(t))j >0 that is given by 

(11) «%a = v 

(12) c\w = -2AVF(a;w), 
with initial condition 

(13) a(0)=0, v(0)=v . 

This process corresponds to a classical particle with Hamiltonian 

if d (a,v) = fef+2AF(a). 

The dynamics of a classical particle in a weak, strongly mixing and homogeneous 
random potential has been studied extensively in the literature, see |14j-[Ti 



We assume that V : R N x Q — > R is measurable such that it is strictly stationary 
in space, i.e., the laws of (V(xi + x ) • • • V{y. n + x )) and (V(xi) • • • V^(x n )) are 
the same for all n G M and x G R . We also assume that 

(14) E(F) = 0, 
and that the realizations of the potential 

(15) V G C 2 (R N ) n W 2 '°°(R N ), P-a.s.. 

Furthermore, we assume that V is strongly mixing in the uniform sense. For 
any r > 0, let C\ and be the a- algebras generated by the random variables 
V(yi), x G .B r := {x G M JV , ||x|| < r} and x G respectively, where is the 
complement of the open ball -B r in K^. We define the uniform mixing coefficient 
between the a-algebras C* and as 

(16) y?(p) = sup{|P(A) -¥(B\A)\, r>0,AeCl., Be C e r+p },p > 0. 
We assume that <p(p) decays faster than any power of p, i.e., for all p > 0, 

(17) sup/Zy^p) < oo. 

p>0 

We introduce the two-point spatial correlation function corresponding to the 
random potential 

(18) R(a) := 4E(F(a)F(0)), 
with Fourier transform -R(v). We assume that 

(19) R G C 00 ^), 
such that 

(20) R does not vanish identically on any H p := {v G R N , v • p = 0}, p G R N . 

In order to observe diffusive behavior of the soliton dynamics for all "macro- 
scopic" times 1/A 2 < Ce\ log/i|/A, as the coupling constant A — > and h — ► 0, we 
need that | \ogh\X — > oo as h — > 0. This is satisfied, for example, if A = m^jot^s , 
for some a G (0, 1). 

We now discuss the limit process as A — > and /i — > 0. Let (v(£))j >0 be 
a diffusion process starting at v(i = 0) = v , where v appears in (jfjj), with 
generator £ whose action on u G Co°(M d \{0}) is given by 

N 

(21) Cu(v) = J2 ^(Ai(v)^.w(v)). 

Here .Dj.,- is the diffusion matrix given by 

l r°° k 

(22) Aj .( k ):=_ / ^^.i?(s — )ds, k G M. . 
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We have the following theorem on the convergence of the momentum of the center 
of mass of the soliton to a Brownian motion on a sphere of constant momentum. 
A stronger result stated for more general nonlinearities, Theorem [5] in Section [2j 
will be proven in Section [5j 

Theorem 2. Consider the NLS equation (TJP in N > 2 dimensions with non- 
linearity f given by (TJj] and potential V satisfying (0)-(O2J). Suppose the initial 
condition satisfies (0|) with ||v || ^ 0, and t/iat there exists a > st/ca t/iat i/ie 
coupling constant A — > as h — > | log /i|A 3/,2+a — ► oo. Then, for any fixed 
and finite T > 0, ine stochastic process (A 2 a(t/A 2 ), v(t/A 2 ))j g ( T ), satisfying (EP 

- nOj) . converges weakly (in law), as \,h — > 0, to (f*v(s)ds, Y.(t))te(o,T): where v 
satisfies l[21\) . 

Actually, in dimensions A^ > 3, we have a stronger result about the convergence 
of the trajectory of the soliton's center of mass to a spatial Brownian motion. 
Let A (x, t, k) satisfy the Liouville equation 

(23) d t (j) x = «%a| 5=Xi v= k ■ V X A + <9tv| 5=x ,v=k ■ V k A . 

with initial condition A (x, 0,k) = o (A 2+/3 x, k), (3 > 0, where 0o is compactly 
supported, twice differentiable in x e IR^, and four times differentiable in k e M^, 
such that its support is contained in the shell 

A(M) := {(x, k) G R 2N , l/M < ||k|| < M}, 

for some M > 1. Let n(x, i, k) be the solution of the spatial diffusion 

(24) d t u = J2dij(\\H)d Xl d x ^ 

with initial condition 

„(x, 0, k) = I o (x, ||k||l)dE(l), 

1 N-l JS"- 1 

where T^-i is the area of a unit sphere S 1 ^ -1 in A^-dimensions, and d£(l) is the 
measure on the unit sphere S N ~ l . Here, the coefficients dij are given by 

dy(||k||) = ^—f ||k||^(||k||l)dE(l), 

1 N-l JS 1 *- 1 

where Xj are the mean-zero solutions of 

N 

<9 fcl (Aj(k)<9 fcj Xj) = -||k||%. 

i,j=l 

We have the following theorem. 

Theorem 3. Consider the NLS equation (TJJ) in N > 3 dimensions with non- 
linearity f given by (TJj] and potential V satisfying Suppose the initial 
condition satisfies wzta ||vo|| 7^ 0, and taat taere exists a > swca taai A — > 
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as h — > wift | log/i|A 1+Q — > oo. T/ien t/iere exists (3 G (0, a/2) swc/i that, for all 
< P < (3, any fixed < t < T < oo, and all compact sets K C A(M), we have 

lim sup |E[</> A (x/A 2+/3 , t/X 2+2fS , k)] - u(x, t, k) | = 0, 

\h->Q (t,x,k)e[t ,T]xK' 

where <f> x satisfies (23\) and u satisfies (24\) . 

The same result holds for more general nonlinearities, see Theorem [6j Section 
[2j We note that the conditions on the convergence of A — > as h — > in Theorem 
[3] are optimal, in the sense that they are the ones needed for the time scale over 
which the soliton behaves like a point particle, 0(\\ogh\/X), is larger than the 
macroscopic time scale over which spatial diffusion is observed, 0(A~ 2_2/3 ), (3 > 0. 

Our analysis of the effective dynamics of the soliton in the external potential 
is based on an extension of [7] to the case where the potential has a random 
character and on the derivation of explicit estimates involving both the spatial 
variation of the external potential h and the coupling constant of the potential 
A; see also [I]- [9]. This analysis relies on three main ingredients. First, using a 
skew-orthogonal (or Lyapunov-Schmidt) decomposition property, we decompose, 
almost surely, the solution of the NLS equation with initial condition close to 
a soliton configuration into a path belonging to the soliton manifold and a part 
describing a fluctuation skew-orthogonal to the manifold. The dynamics on the 
soliton manifold is obtained by the skew-orthogonal projection of the Hamiltonian 
flow generated by the NLS equation in a small tubular neighbourhood of the 
soliton manifold onto the latter. As for the fluctuation, we control its if ^norm 
almost surely using an approximate Lyapunov functional. To study the limiting 
dynamics in the weak-coupling/space-adiabatic limit, we rely on explicit control 
of the difference between the effective dynamics of the soliton and a classical 
particle in the random potential and on the results of [H],[l5] and [17] on the 
stochastic acceleration of a classical particle in a random potential. 

The organization of this paper is as follows. In Section [2J we list our general as- 
sumptions and state our main results for general nonlinearities: Theorem [Hon the 
long-time dynamics of the soliton in the random potential, and Theorems [5] and 
[6] on the limiting diffusive dynamics in the weak-coupling/space-adiabatic limit. 
In Section [31 we recall basic properties of the NLS equation, which we will use in 
the analysis. We then prove Theorem H] in Section HI In Section [5] we prove The- 
orems [5] and [6] by showing first that the limiting effective dynamics of the soliton 
converges almost surely to that of a classical particle in a random potential and 
then applying the results of [2] , [T5] and [T7] on the motion of a classical particle 
governed by a weakly random Hamiltonian flow. For sake of completeness, we 
discuss the well-posedness of the generalized NLS with random interactions and 
potential in Appendix A. We also prove some technical statements in Appendix 
B. 

1.3. Notation. 



• Given x G M. N , we denote x = ||x|| := y^ i=1 i t 2 and x := ^j. For 

x, y G R^, we denote their inner product by x • y := J2f =1 Xiyi. 

• L P (I) denotes the standard Lebesgue space, 1 < p < oo, with norm 

\\f\\i, = (Jdx\f(x)\*)*, fetf(I),p<oo, 
H/ll^ = ess -sup(|/|), /GL°°(7). 
We also define 

\\f\\LP(I,Li(J)) ■■= || \\lp(I)- 

For 1 < p < oo, p' is the conjugate of p, i.e., l/p+ l/p' = 1. We denote by 
(•, •) the scalar product in L 2 (IR Ar ), 

(u,u) = Re / uv, it, f G L 2 (M 7V ), 

as well as its extension by duality to Y x Y', where Y and Y' are complete 
metric spaces such that Y L 2 Y', with dense embedding 

• Given the multi-index a = (ai, ■ ■ ■ ,«jv) G N^, we denote |a| = YliLi a i- 
Furthermore, d" := d"l • ■ ■ d%£. 

• For 1 < p < oo and s G N, the (complex) Sobolev space is given by 

W S ' P (R N ) := {u G S'{R N ) : <9> G ^(R^), |a| < s}, 

where S'(M. N ) is the space of tempered distributions. We equip W s ' p with 
the norm 

\\u\\ w ., P = W d x u hp, 

a,\a\<s 

which makes it a Banach space. Moreover, W~ s,p ' is the dual of W s ' p . We 
use the shorthand W s ' 2 = H s . 

• Given / and g real functions on R N , we denote their convolution by *, 

/ * 0(x) := y dy /(x - y)#(y). 

• We say that a real function / G C 6 1,1,2 (R JV x [0, +oo) x M w \{0}) if it 
satisfies 

sup sup |c^9f3£/(x, £, k)| < oo. 

N<1,/3<1,|7|<2 (x,t,k)6]R iv x[0,+oo)xIR jv \{0} 

2. General nonlinearities 

The model. We now list our general assumptions and discuss models where 
are satisfied. 



(Al) Nonlinearity. The nonlinearity f — fi + • • ■ + fk such that each 

j 3 G C 2 (H\R N ; C), H~\R N ; C)), j = 1, • • • , k. 

We assume that there exists 

2N 

r i e &W^2 l ([2 '°° ] if iV = 1) ' 
such that VM > 0, 3 a finite constant Cj(M) such that 

\\fj( u ) ~ fji^Wjf'j ^ Cj(M)\\u — v\\ L r j: 

Vu,v G H\R N ;C), \\u\\ H i + \\v\\ H i < M. Furthermore, 

hafj(u)u = 

almost everywhere on R N ,\/u G i/ 1 (]R Ar ; C). 

Let F := Ylj=i Fj: "where Fj satisfies fj = Fj, and the prime stands for 
the Frechet derivative. We have that 

3Fj G (^(/^(R^C),^). 

For every M > 0, there exists a positive constant C(M) and e G (0,1) 
such that 

F(w) < ^^Ikll^ + C(M), Vu G tf 1 ^) 
such that ||m|| L 2 < M. Furthermore, 

SUp ||F"(u)|| B (Hi 1 H-i) < 00 

\\u\\ Hl <M 

sup ||F'"(u)|| H i_ B ( H i )H -i) < 00, 

\\u\\ Hl <M 

where B denotes the space of bounded operators. 
(A2) Symmetries. The nonlinearity F satisfies F(T-) = F(-), where T is a 
translation 

Tf : u(x) -> u(x-a), a G R N , 

or a rotation 

T£ : u(x) -> u(i2 _1 x), i? G SO(N), 
or a gauge transformation 

T» : u(x) -> e i T«(x), 7 e [0,2tt), 

or a boost 

T h v : «(x) -> et v ' x M (x), vGl N 
or a complex conjugation 

T c : u(x) -> u(x). 
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(A3) Solitary Wave. ]/cl such that V/i G /, the nonlinear eigenvalue prob- 
lem 

(25) (-A + /i)^-/(^) = 

has a positive, spherically symmetric solution r/^ G L 2 (R N ) fl C 2 (WL N ), 
such that 

Hllxll^VlU 2 + IIII x II 2 IIv?7 m ||||l2 + iiiixii 2 ^^!!^ < oo,V/i g /. 

(A4) Orbital Stability. The solution appearing in assumption (A4) satisfies 

d M m(» > 0, V/i G /, 

where m(/i) := | J <ix r^ 2 is the "mass" of the soliton. 
(A5) Null Space Condition. We define 

C„ := -A + fi — f(v»), 

which is the Frechet derivative of the map ip — > (—A+fi)i/j — f(ip) eval- 
uated at r/^. For all [A G I, the null space 

■Af(A») = spaniirj^ d^rj^j = 1, • • • , JV}. 

(Bl) External Random Potential. The external potential is a real measurable 
function on x K x -> R such that 

V^(x,t;a;) =^(/ix,*;w), he (0,1], 

and it realizations 

V(x, t; uo) G VF 1,0O (R, C 2 ^) n W 2 ' co {R N )), P - a.s.. 

(B2) Time-Independent, Homogeneous and Strongly Mixing Random Potential. 
The external potential is time- independent, 

V(a,t;u) = V(a;uj), 

and V is a homogeneous random potential that is strongly mixing in the 
uniform sense, satisfying conditions (I7j)- (|2"U1) . Subsection 11.21 

Remark 1. Assumptions (Al) and (Bl) are sufficient to establish almost surely 
the global well-posedness of the NLS equation (E]) in H 1 (see Appendix A). More- 
over, (Al) implies 

\F{u + v)- F{u) - (F'(u),v)\ < C(M)\\v\\ 2 H i 

\F(u + v) - F(u) - (F'(u),v) - ±(F»(u)v,v)\ < C(M)\\vf H1 

\\F'(u + v)- F'{u) - F"(u)v\\ H -i < C{M)\\vf m , 

for any u G H 2 (R N ) and v G H 1 ^) such that \\u\\ H i + \\v\\ H i < M. 
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Remark 2. We now discuss typical nonlinearities where the various assumptions 
are satisfied. Assumptions (Al) and (A2) are satisfied, for example, if 

F(u) = \J d * G iH 2 ) + w * M 2 , 

where G(r) = JJ" dsg(s), such that g G C 2 (R + ) with g(s) < C(l + s a ), a G [0, |), 
19*0001 <C(l + S "- fc ),A; = 0,l,2,gG [0, if N > 3, q G [O,oo), if N = 

1,2. Furthermore, W G L p + L°°,p > f ,p > 1, such that max(0, G L r (M^) + 
L°°(R Ar ),r > y, > 1, if N > 2, and W is spherically symmetric, see [27]. 
Furthermore, (A3) is satisfied for local nonlinearities if 

— oo < limgfs) < )U 
-oo < lim s' a g(s) < C, 
where < a < 2/(N - 2), when N > 2 and a G (0, oo) if N — 1, 2, such that 

r< 

3( > 0, with / dsg(s) > 
Jo 

see [211 (22J [27J . Assumption (A3) is also satisfied for nonlocal nonlinearities if, 
in addition to the above, 

N 

W G L\ oc , q > —, and W — > as ||x|| — > oo; 

see J27J [201 [21 [3] . Assumption (A4) implies orbital stability of the solitary wave 
solution, |23j . /£ zs satisfied, in particular, for local nonlinearities 

/(VO = IV-ftM < ^- 

Assumption (A5) is satisfied for local nonlinearities if 

g\s)+g"{s)s 2 >^ 

orifN=l, 011231. 

2.2. Statement of the main theorems for general nonlinearities. We now 

state the main results of this paper. 

Theorem 4. Consider the NLS equation (QJ) with initial condition satisfying |6p. 
Suppose that (Al)-(A5) and (Bl) hold. Then the result of Theorem^ hold. 

We now discuss convergence of the dynamics of the center of mass of the soliton 
to Brownian motion in the weak-coupling/space-adiabatic limit. 

Theorem 5. Consider the NLS equation (QJ) in dimensions N > 2 and with 
initial condition satisfying (0|) and \\vq\\ ^ 0. Assume hypotheses (A1)-(A5), 
(Bl) and (B2) hold. Suppose in addition that there exists a > such that A — > 
as h — > with \ log/i|A 3 / 2+Q — > oo. Then the result of Theorem^ hold. 
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We also have the following stronger result for the convergence of the trajectory 
of the soliton to spatial Brownian motion in the weak-coupling/space-adiabatic 
limit in dimensions N > 3. 

Theorem 6. Consider the NLS equation (T7J) in dimensions N > 3 and with 
initial condition satisfying and ||vo|| 7^ 0. Assume hypotheses (A1)-(A5), 
(Bl) and (B2) hold. Suppose in addition that there exists a > such that A — > 
as h — > with | log/i|A 1+Q — > 00. T/ien i/ie result of Theorem^ hold. 

We prove Theorem H] in Section HJ and Theorems [5] and [6] in Section [5j Before 
doing so, we recall some useful properties of the NLS equation. 

3. Mathematical preliminary 

In this section, we recall some basic properties of the nonlinear Schrodinger 
equation ([T]), see for example [I]. 

3.1. Hamiltonian structure. The phase space for the NLS equation (pQ) is cho- 
sen to be H l (M. N ; C). The space if 1 (IR iV ; C) has a real inner product (Riemannian 
metric) 

(26) (u,v) := Re / dx uv 



for u,v G if 1 (IR iV ;C). The tangent space at an element ip G H 1 is T^H 1 = H 1 . 
On i? 1 (R JV ; C), one can define a symplectic 2-form 

(27) E(u,v) := Im ^ dxuv= (u,iv). 

The Hamiltonian functional corresponding to the NLS equation (CQ) is 

(28) H X W) := i j dx I V^| 2 + ±jv h \4,\ 2 - F(V»). 



Using the correspondence 

C) < ► /^(R^; R) © iJ 1 ) 

■0 « — > (Re?/> Inr0) 



)JV. 



-1 



J. 



where J := f ^ J J is the complex structure on if^R^; R 2 ), the NLS equation 
can be written as 

W = JH' X ^). 

We note that since the Hamiltonian functional H\ defined in ( 128]) is nonau- 
tonomous, energy is not conserved. For ip G H 1 satisfying ([!]), we have that 



d t H x {^) = ^ j dx(d t V h M 2 , P 



a.s., 
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see Appendix B for a proof of this statement. Yet, H\ is almost surely invariant 
under global gauge transformations, 

H x (e^) =H X {$), P-a.s., 

and the associated conserved quantity is the "charge" 

Nty) -=\j*x. |^| 2 . 

The assumption d M m(/x) > 0, where m(/i) is defined in (A4), implies that 77^ 
appearing in assumption (A3) is a local minimizer of H\ =0 (tp) restricted to the 
balls B m := {ip G H 1 : N(ip) = m}, for m > 0; see [23]. They are critical points 
of the action functional 

(29) £ M (y,) := \ J dx (|V^| 2 + /x|^| 2 ) - Fty), 

where /x = /i(m) is a Lagrange multiplier. 

3.2. Soliton Manifold. Let 77^ be an (orbitally stable) soliton solution of (121)]) . 
and define the soliton manifold as 

M s := K := T av7 r/^, a = (a,v, 7 ,/i) e^xl^x [0, 2vr) x /}, 

where I appears in assumption (A3), and the combined transformation T av7 is 
given by 

V^av 7 := T av ^ = e^v-(x-a) +7 )^ (x _ a)> 

where v, a G M. N and 7 G [0, 2tt). Note that if the nonlinearity / satisfies /'(0) = 0, 
then J C M + . 

The tangent space to at 77^ G A^ s is given by 

T^M S = span{E t , E g , E b , E s }, 

where 



E t 
E 9 



V a Tf 77 M | a=0 = -V77, 
2V v T^| v =o = ^x?/V 



F s := 9 M 7/ M . 

In the following, we denote by 

ej := -d Xj , j = ,N, 
e j+N :=ixj, j — 1, ■ ■ ■ , N, 

627V+1 := h 

(30) e 2N +2 ■= <9 M , 

which, when acting on 77^ G A^s, generate the basis vectors {e a ri a }l t I l~l 2 of T Vcr M. s . 



(31) 
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The soliton manifold M. s inherits a symplectic structure from (H 1 , S). For 
a = (a, v, 7, fi) eR N xR N x [0, 2vr) x 7, the matrices 

- = P 7 _1 P 

where P a is the L 2 -orthogonal projection onto T^Ais, define the induced sym- 
plectic structure onX s . Explicitly, we have 

^•a\T v<T M a '■ = {{ e af](T^ e l3' r la)}l<a,l3<2N+2 

( -m(fi)l NxN -lvm'(jj)\ 

m(fx)l NxN am'(/i) 

ra'O) 

V |v T m / (/i) -a T m'(yu) -m'(/i) 7 

where Inxn is the N x N identity matrix, and (-) T stands for the transpose of a 
vector in Mr. One can easily show that if 9 At m(/x) > 0, then S CT is invertible. 

3.3. Group structure. The anti-selfadjoint operators {e a } a= i ) ... ,27V+i defined 
in fl30l form the generators of the Lie algebra g corresponding to the Heisenberg 
group H 2N+1 , where the latter is given by 

(a, v,7) • (a , v ,7 ) = (a , v ,7 ), 

with a" = a + a', v" = v + v', and 7" = 7' + 7 + |v • a'. Elements of g satisfy the 
commutation relations 

(32) [e», e j+N ] = -e 2 N+Aj, i,j = 1,- - ,N, 
and the rest of the commutators are zero. 

3.4. Zero modes. The solitary wave solutions transform covariantly under trans- 
lations and gauge transformations, i.e., 

^(T'JT^) = 

for all a e M N and 7 e [0, 2vr). Here, the prime stands for the Frechet derivative. 
There are zero modes of the Hessian, 

(33) C,j, :— —A + /x — /'(%), 
associated to these symmetries. One can show that 

iCp : T^M, -> 2^.M S 

with {iC^) 2 X = 0, for any vector X e T Vfi M s . 

To see this, differentiate £' (T^T]^) = with respect to a and set a to zero, 
which gives 

(34) £"(^)V a ^(x _ a )| a=0 = CpEf = 0. 

Similarly, differentiating £' (T-frjn) = with respect to 7 and setting 7 to zero 
gives 

(35) C^Eg = 0. 
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Using (1251) . we have 

(36) L^E h = iE t , 
and 

(37) C^E S = iE g . 

3.5. Skew-Orthogonal Decomposition. Consider the manifold A4' s = {7/0-, a G 
E }, E = R N xR N x [0, 2tt) x J , where I C J\<9/ is bounded. We define the 
5 neighbourhood of A4' s in i? 1 as 

(38) U s := {ip G #\ inf ||V - v*\\w < 5}. 

Then, for 5 small enough and for all ip G t/j, there exists a unique <t(V>) G 
C^t/^E) such that 

- = (ip- r] aW ,iX) = 0, 

for all X G T v M. s . For a proof of this statement, we refer the reader to [3]. 

Remark 3. TTie group element T av7 G H 27V+1 zs ^wen fry 

T av7 = e'^e^e^. 

It follows from UM) that T^LYT av7 G g i/ K G g. Furthermore, it follows from 
translational invariance that S(T av7 w, T av7 f ) = E(u,v), foru,v G L 2 . Therefore, 

Z(w,Y Va )=Z(w',Y'ri a/ ) = 0, 

VY G g, where Y' = T~^YT mi e g, w' = T av \w, and = T av \r] a . 

4. Long-time dynamics in a random potential 

In this section, we study the long-time dynamics of the center of mass of the 
soliton, and we prove Theorem HJ The main ingredient of our analysis is an 
extension of the analysis of solitary wave dynamics in time-dependent potentials, 
[7], to the case of random potentials; see also [3J-[6], [8] and [9] for the time- 
independent case. Furthermore, we improve slightly the error bounds in [7j by 
obtaining estimates that depend both on the coupling constant A and the scale 
of the spatial variation of the external potential, h. 

The weak solution of ([T]) with initial condition cf) is given by the Duhamel 
formula 

(39) i){t) = U(t,Q)(j)-i\ f dsU(t,s)V h (x,t)i/j(s)+i [ dsU(t,s)f(s,ip(s)), 

Jo Jo 

where U(t, 0) = e lAt is the unitary operator corresponding to free time evolution. 
One can show that it is unique, and that it is P-a.s. in H 1 ; see Appendix A 
for a proof. It follows from (139]) that if) is cj-measurable. We have the following 
proposition. 
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Proposition 1. Consider (T7J) with initial condition satisfying (E|). Suppose as- 
sumptions (A1)-(A5) and (Bl) hold. Then there exists h > and positive 
constants C and C, such that, for all h G (0,/iq), any fixed e G (0,1), and all 
times t G [0, Ce\ log/i|/A/i), we have that 

V>(x,*) = 7fc(t)( x ) + w(x,t), LJEtl, 

with 

sup sup IIHIi?" 1 < C/i 1_ 2 ; 

weH te[0,Ce|logft|/Aft.) 

uniformly in A G (/i 1 " 6 , 1], an<i t/ie parameters a, v, 7 and /i satisfy the effective 
differential equations, for uefl and t G [0, Ce| log/i|/A/i), 

(40) <9 t a = v + 0(/i 2 ~ £ ) 

(41) d t v = -2AVy„(a, t) + 0(/i 2 - e ) 

(42) da = ii- AV^ (a, t) + ^ 2 + 0(/i 2 - £ ) 

(43) 9 4 /i = 0(/i 2 - e ), 

™^ ||a(0) - a ||, ||v(0) - v ||, |t(0) - 7 |, |//(0) - /i \ = 0(h). 

The proof of Proposition [1] is given in Subsection 14.31 We first find the 
reparametrized equations of motion corresponding to the restriction of the Hamil- 
tonian flow generated by the NLS equation to the soliton manifold. 

Remark 4. In the special case of local nonlinearities in dimension N = 1, and 
under the strong assumption on the initial condition, 

tjj(t = 0) = ri ao , 

one obtains a slightly better estimate on the H 1 norm of the fluctuation w in 
Proposition [21 

sup sup HHIh 1 — Ch 2 ~i, 

uieTi te[0,Ce\logh\/Xh) 

see, for example, |9J for a discussion of the cubic one dimensional NLS equation 
with a slowly varying time-independent potential. Here, we choose to work with 
general nonlinearities in all dimensions under weaker hypotheses for the initial 
condition. 

4.1. Reparametrized equations of motion. Suppose ip satisfies the initial 
value problem (TTJ such that, P-almost surely, if) G Us C H 1 , where Us is defined 
in (1551) . Subsection 13.51 By the skew-orthogonal decomposition, there exists a 
unique 

a = a(ip) = (a, v, 7, //) G S = R N x R N x [0, 2tt) x I 
and w' G H 1 such that 

lf> = T] a + w', \/0J G Q, 
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with 

w' JL J-%„M S . 

For uj E Q, let 

(44) m := T a "!^ = ^ + w, 

where w = T~^ w'. 

We introduce the coefficients 

Cj := dtdj - vj, c N+j := --d t Vj - A9 a . J .V r h (a, t), j = 1, • • • , N, 
1 i 

(45) c 2 n+i := fi- -v 2 + -d t a ■ v - XV h (a, t) - d t j, c 2 n+2 ■= 
We denote by 

(46) |c| := sup |cj|, 

ie{l,-,2iV+2} 

and 

C(|c|, w, h) := sup{|c| HwIIh 1 + Xh 2 + ||if||^i}. 
We have the following lemma. 

Lemma 1. Consider the NLS equation (QJ) with initial condition Suppose as- 
sumptions (A1)-(A5) and (Bl) hold, and that for uo G Q, ip{t) G Ug fort G [0,T], 
where 5 and T are independent of uj G Q. Then the parameter a = (a, v, 7, /1), as 
given above, satisfy 

d t dj = vj + 0(C(\c\,w, h)) 

d tVj = -2d Xj XV h {a,t;u)+0{C{\c\,w,h)) 

da = n-- ||v|| 2 + -dtSL ■ v - XV h (a, t; uj) + 0(C(\c\, w, h)) 
d t tM = 0(C(\c\,w,h)), 
forte(0,T),j = l,---N. 

Proof. We first find the equations of motion in the center of mass reference frame. 
For uj G £1, we differentiate 

u(x, t) = T a ->(x, t) = e-^ v ' x ^V(x + a) 
with respect to t, and use the fact that 

llvll 2 

e -| (vx+7) A ^( x + a ) = A«(x, t)+iv- V«(x, t) - 



e -^^/(^(x+a,t)) = /( M (x,t)). 
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We have that 

2iV+l 

(47) d t u = — A + fi)u — /(«)) + c j e i M — ^v u 

3=1 

where 

7£y := AVfe(x + a, t; u) - A\4(a, t; u) - AV\4(a, t; u) ■ x, 

and ej and Cj are as defined in fl30l) and f|45l) respectively. Equivalently, we have 
that 

2N+1 

d t u = —i£'^{u) + c i e j M — i'R-v u , 

3=1 

where £ M appears in (1251) . 

We now use f|4T|) and the skew-orthogonal decomposition to find equations for 
the parameters a = (a, v, 7, /i) and w, for w6fi. Since 

we have that 

where £ M = (-A+/x-/' (^)) = ^'(^) and iV» = f^+w)-f^)-f'^)w. 
Substituting this back in (T4TI) gives 

2AT+1 2iV+2 

d t w = (—iCp + ^2 c J e 3 ~ iR-v)w + N^w) + - iRyf]^ 

j=i 3=1 

Now, we know that (X,iw) = for all X e T V A4 S , uj E £1, see Remark [31 It 
follows that 

<9 t (X, wu) = <9 4 /i (d^X, iw) + (X, id t w) = 0. 
Therefore, for u G O, we have that 

2iV+l 

= -(iX,iCfj,w) - (iX,iTZ v (r)f, + w)) + (iX,N^(w)) + (iX, c^-to) 

2iV+2 

i=i 

It follows from CT-CT that (iXJC^w) = 0. Together with 
[ ei ,i]=0, e* = -e v j = !,■■■ ,2iV + 2, 

we have that 

27V+2 2W+2 

(48) C ^ X ' iw ) = ~^ X ' + ^) + ^H) + J2 C ^ X ' z ' e i^)' 
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for uo G £1. Now, choosing X = Ek, where Ek, k G {1, • • • , 2N + 2}, is a basis 
vector of T^A^ gives 

2N+2 2N+2 

^2 ( 5 )fci c j = ( E k, N^w) + K v (w + rjf,)} + ^ Cj(iejE k , w), u G O, 

3=1 3=1 



where 5 appears in ( 13T1) . Replacing the definition of Cj, appearing in ( 1451) . in 
181) . and using the fact that 



(x,-^, i^y^) = 0, ri^x) = r/ M (||x||), (r^, zfty^) = 0, 
gives, for uj G Q, 

. 2N+2 

(49) <9 t a fc = w fc + ^J^(( ix kV^ N^w) - %R v w) + c^efi^, w)), 
(50) 

2 2iV+2 

9 t u fe = -2d Xk XV h (a.,t;uj)-\ -—{{d Xk r]^,N^(w)+TZ v w)- c j (ie j d Xk r] fl ,w) + (d Xk r]^,'JZ v r]^)), 

r7 HA t J . =1 

% =/i - ^||v|| 2 + ~<9 t a- v - AV^a, t; u;) - ^ ((<Vfc, -^(w) + 

27V+2 

(51) - ^ Cjiiejd^w)), 

3=1 

- 2iV+2 

(52) <9 t /i = —^.(iTj^N^w) -iK v w) - ^ Cjie^w). 

The claim of the lemma follows directly from Assumptions (Al) and (A3), 
which imply that 

sup HT^y-EfcHi^ = 0(\h 2 ), 



and 



sup ||iV M (w) ll^i < C sup|H|^i 5 



for sup^gQ HwIIh 1 < 1 an d some constant C that is independent of /i and A. 

□ 



20 

4.2. Control of the fluctuation. We use an approximate Lyapunov functional 
to obtain an explicit control on sup we ^ an d su P^en l c l- This approach dates 

back to [24], and has been used, with various generalizations, in [I] -IS]. We define 
the Lyapunov functional 

(53) C^u,v) :=£»-£», u,veH\R N ), 

where £ M is given in (1291) . We proceed by estimating upper and lower bounds for 
Cfj,(u, rjfj), uefl, where u = rj^ + w is defined in (144"I) . 

Lemma 2. Suppose the hypotheses of Lemma[J\hold. Then there exists a constant 
C independent of h such that 

(54) sup IdtC^u,^)! < C sup (Xh 2 \\w\\ H i + (|c| + Xh + \\w\\ 2 m )\\w\\ 2 H1 ) , 

where \c\ appears in $fi\) , uniformly in t G (0, T) and A G (/i 1_e , 1]. 

Proof. For a; G fi, we have 
1 



|2 



= - / dx |Vm| 2 + //[up - F(u) 

= ^A(T a -» + \^\\T av \nh - \ Jdx V h \T- v \iP\'< 
By translational symmetry, we have, for ui G Q, 

h\\» = W\\h, J d*v h \u\ 2 = j wi 2 , 

where V^°(x) := V h (x + x ), x G R N . Moreover, for ueO, 

^A(T a -» = ^A(^) + ^(i|| V || 2 +>u)||^||| 2 -i V -<^,V^) + ^ J dx(V-*-V h M 



and hence 



(55) ^( M ) = ^ A W + i(i||v|| 2 + ^)||^||i 2 -iv(#,V^)-^ jdxV h \^\ 2 . 

We have the following relationships for the rate of change of field energy and 
momenta, whose proof we give in Appendix B. For uj G Q, 

(56) d t ^J dx v h \^\ 2 = (VW,V^) + i<V,$W), 

(57) ^tf A (^) = ^,$W) 
and Ehrenfest's theorem 

(58) ^(#,V^) = -(^,AV^). 

One can use (JT]) to formally derive (|56j) . (!5T|) and (I581) . Rigorously, one can prove 
them using a regularization scheme and a limiting procedure, and we refer the 
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reader to the Appendix B for a proof of these statements. Furthermore, it follows 
from gauge invariance of ([T]) that the charge is P-a.s. conserved, 

(59) $|MU» = 0, cuEQ, 

(Proposition [2] in Appendix A) . Differentiating fl55|) with respect to t and using 
gives, for to G Q, 



"1 s~) 1 1 

d t £,{u) = dtHxty) + - (-^ + fyi) Ml - -d t v • W> + -v • (V, AVW> 

= 2 (^y^ + MINI!* - ■ v ^ + 2 v ' XVVhij) ~ < AVV ^' V ^ 



(60) 



- (~<*($v + 2AW,>, Vu), 



where we have used w(x, i) = e^2 v ' x+ T)^(x + a, t) and translation invariance of 
the integral in the last line. Furthermore, it follows from (|29|) that 

which, together with fl53|) and fl60l) . implies 

(61) ifo) = ^MINIi* - |foX) - {i{\d t v + AW,>,Vu>, 

for a; G f2. We now estimate the right-hand side of (16TT) . Since (zX, w) = for all 
X G T^M. S and co> G f2, we have that 

Hi* - IMlia = \\ w \\h, 

and hence 

(62) ^MIHIi»- Wi») = o(|c|IHIi>)- 

To estimate the second term in the right-hand side of flBTl) . we replace u = 
77^ + w, and use the fact that (iE g , w) = (iV^, w) = 0, and (igrj^, Vt^) = for 
all real # G L°°(lR iV ) and u G H. We have 

(^tv+AV^>, V«) = (iAW>, V^> + (zAW h a r/ M , Vu>)+<i(^v+AW fc a )w, V»), 

for w£(l. Adding and subtracting the quantity 

AV14(a, t; oj) ■ (iw, Vr]„) = XVV h {a, t; u) ■ Vw) = 
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gives 

(z(~<9 t v + AW,*)u, Vu) = (i^tv + A W^(a, *; Vw) + ((AW h a - AW fe (a, *; w))iu;, Vio) 

+ ((AVV£ - AW fc (a, t; u))i Vlt , Vw) + ((AW? - AVV^a, *; w))iu;, V^), 

for u; E fi. It follows from Lemma [T] that the first term of the above equation is 
of order 0(|c|||w||^-i + A/i 2 ||w||^i + while Assumption (Bl) implies that 

the second term is of order 0(A/i||w||^i). Assumptions (Bl) and (A3) imply that 
the third and forth terms are of order 0(A/i 2 ||u>||#i). 

□ 

The object of the next lemma is to provide a lower bound for sup wg ^ \C^(u, rj^) | . 
Let 

X M := {w E H\R N ) : (w,iX) = 0, WX E T^M S ). 
It follows from the coercivity property of that there exists a positive constant 

(63) p := inf (w,£uw) > 0, 
(see Appendix D in [4] for a proof of this statement). 

Lemma 3. Suppose the hypotheses of LemmaUl hold. Then there exists positive 
constants p and C independent of h and uo such that, for sup wg Q < 1 and 

uniformly in t G [0, T], 

(64) sup|C M (u,7/ M )| > sup^HHItfi -C\H\%i], 
where p is defined in / Tffgj) . 

Proof. For uo G Q, we expand S^u) around rj^, which is a critical point of £ M . 

(65) 8^ + w) = SM + ±(w, Z» + Rf{w), 
where 

Rf{w) = F(r )tl + w)-F{r ] ^ - (F'(^),w) - ±(F"{ Vll )w,w). 

It follows from Assumption (Al) that 

sup \Rf\w)\ < Csup 

for sup^gQ ||w||^i < 1, where C > is independent of t G R. Moreover, the 
coercivity property (1631) implies 



SUp(u>, CfjW) > pSUp 

and hence 

sup IC^tt,?^)! = supj^w) -8^)1 > sup{^p||w||^i 
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for sup w6iJ \\w\\ H i < 1. □ 

We now use the upper and lower bounds on the Lyapunov functional to obtain 
an upper bound on sup wg Q 

Lemma 4. Suppose (A1)-(A5) and (Bl) hold. Let ip satisfy (Tip and, foru G tt, 
n, 77^,-0; be given as above. For h -C 1, choose T G IR + suc/i £/ta£, /or w G Q, 
ip(t) G L^, t G [0, T], where Us is defined in / fffgj) . Subsection \3.5[ Fix e G (0, 1), 
and choose t G [0,T] stzc/i taai sup^g^ ||w(to) < h 2 ~ e . Then, for h small 
enough, there exist absolute constants C\ > 1 and C2 > 0, which are independent 
of h and e, such that 

sup sup < Ci(^ 2 + sup ||w(to)||^i)) 

wer2*e[t ,to+C"2/A/i] w ea 

sup sup |c(t)| < Ci(h 2 + sup ||w(io)||lfi), 

a>efi*e[io,io+C2/A/i] weQ 

uniformly in A G {h l ~ e ) 1]. 

Proof. It follows from Lemma [2] that, for t >to 

sup |C M (^ + u>(t),n M )| < sup{|C A1 (n At + w(t ),n M )| 

+ C(t - to)(|c||k(t)||^ + A/i 2 ||«;(t)||^ + A/i||n;(t)||^ + |MI^)}- 

Expanding £^(77^ +w(to)) around n M and using Assumption (Al) gives the upper 
bound 

|C M (n M + w(t ), 7)^)1 < C\\w(t )\\ 2 H i, for \\w(t )\\ H i < 1 and w G O, 
where C is a constant independent of /i, e and A. Therefore, 

sup|C M (7fc+w(t),7fc)| < C sup ||w(to)||^i+C(t-t ) sup(A/i 2 ||w(t)|| H i + (|c|+A/i+||w||^i)||u;(t) \\%i), 

u£f2 w£f2 wSf2 

for some constant C independent of h, e and A. Together with LemmalU it follows 
that 

-supp\\w(t)\\ 2 H i < C sup \\w{t ) \\ 2 m +C(t-t ) sup(A/i 2 ||n;(t) \\ H i + (|c| + Xh + \\w\\ 2 H x)\\w(t)\\ 2 m ) 
+ Csuv\\w(t)\\ 3 H i 

where p appears in floHj) . Equivalently, there exists a positive constant C inde- 
pendent of h, e and A, such that 

C sup\\w(t)\\ 2 H i < sup \\w{t ) \\ 2 H i + (t - t ) sup(\h 2 \\w(t)\\ H i + (|c| + \h+ |MlfrOIK*)il|i) 
+ sup 
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For 



C 

2(A/i + |c| + sup^ en H^ll^i) 
Csup \\w(t)\\m < sup{|k(t )|||i + ^h\\w(t)\\ m + ^\\w(t)\\ 2 Hl + \\w(t)f Hl }. 



Using the fact that 



we have 



h\\wmm<lh a + ^\\w(t)\\m, 



C C 

sup ||w(io)||#i + ~jh 2 ~ "T SU P imOllffi + SU P ll^Wlllf 1 > 0- 

Let ?/o := sup wgT! ||w(t )||fl-i, 2/ := sup wg?! sup 4e[t0jto+T] \\w{t)\\ H i, and /(j/) = - 
j?/ 2 + Vo + jh 2 . For /i<l, the function intersects the x-axis in a point such 
that ^ < ci(h 2 + jjq), where c\ is a positive constant independent of h and e. It 
follows, if ?/o < 2/*, that y < y* for t G [to,*o + r ]- Substituting back in (149H52H 
and using (1431) and (HoT) gives 

sup|c| < C2O 2 + 2/o)> 
wen 

for some positive constant C 2 that is independent of h, e and A. It follows that for 
h small enough, there exists positive constants C\ and C2 which are independent 
of h and e, such that 



sup sup \\w(t)\\ 2 Hl < Ci(h 2 + sup\\w(t )\\ 2 H i) 

ueTi te[to,t +^] uen 

sup sup \c(t)\ < C\(h 2 + sup \\w(t )\\ 2 H i)i 

uniformly in u; G f2 and A G (/i 1_e , 1]. □ 

4.3. Proof of Proposition Q], We are now in a position to prove Proposition 
[T]by iterating Lemma H on time intervals of order O(tt) and using the result of 
Lemma [TJ 

Proof of Proposition^ Fix e G (0,1). For u G and /i small enough, let T* 
be the maximal time for which the skew-orthogonal decomposition is possible. 
Consider the interval 

[0, T] = [t , h] U [h, t 2 ) U • • • U [t n _ x , t n ] C [0, T% 

such that 

C 2 

= t < h < ■ ■ ■ < t n = T, (t i+ i -U) < — , i = 0, • • • , n - 1, 

An, 



25 



where Ci appears is Lemma HI We will choose n£N depending on h and e later. 
Let 

\c\i := sup sup \c(t)\, 

ueU te[ti,t i+1 ] 

Hi := sup sup ||«;(£)||#i, i = 0, • • • , n — 1. 

wen t6[ti,tj+i] 

Note that yo < h and |c|o < C7t, for some constant C independent of /i, A and e. 
Iterating the application of Lemma H] we have 



|c|„ < CC^ 1 ^ 2 . 

We now choose n such that C" +1 /i 2 < /i 2_e . This implies 

log h 



n + 1 < -e 



logC; 



Therefore, for i e [0, e ° 2 |logfc| 



logCi Ah J ' 

2-e 



sup |K*)||^1 < hf 

sup|c(t)| < C/i 2 " e . 

u>eTi 

The effective equations for the parameters on the soliton manifold follow from 
the above estimates and Lemma [TJ Furthermore, it follows from ([6]) and the 
skew-orthogonal property that 

||a(0) - ao||, ||v(0) - v ||, | 7 (0) - 7 o|, HO) - /i | = 0(h). 

□ 



Equipped with the above results, we are in position to prove Theorem HI 



4.4. Proof of Theorem HI In Proposition [T], the time-dependent parameters 
a(t) = (a(t), v(t), j(t), fJ>(t)) of the soliton solution satisfying equations (JlT?l) - (jl5|) 
are defined for uj G £1 (a dense set of Q). We extend them to be random variables 
of all the realization space Q. 
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Proof of Theorem^ We extend (a(t), v(t), j{t), ^{t)) <t<Ce\iogh\/xh appearing in 
Proposition [I] to random variables on Q by requiring that they satisfy 

^ 2N+2 

d t a k =v k + — ({ixk-q^N^w) - iK v w) + ^ c j( e j x kVit>w)), 

j=i 

2 

d t v k = - 2d Xk \V h (&,t;u) + jj^{{d Xk r]^ N^(w) +1Z v w) 

2N+2 

dtl =H - ^||v|| 2 + -d t a ■ v - XV h (a, t; u) - f H^T^yU 9 ** 7 ^ N v( w ) + n v{w + %)) 

2N+2 

. 2V+2 

where 

with the same intial conditions appearing in Proposition [H which are determined 
by the skew-orthogonal decomposition property. Since V is also cu-measurable, 
o~(t) = (a(t), v(t), 7(£), (J>(t)) and r) a (t) are cu-measurable for <t< Ce\logh\/\h. 
Note that we have from Proposition [1] that, for < t < Ce\ \ogh\/Xh, 

nm)-^(t)\\m<ch 2 -\ 

which completes the proof of the theorem. □ 

5. A LIMIT THEOREM IN THE WEAK-COUPLINC/SPACE-ADIABATIC REGIME 

In this subsection, we prove Theorems [5] and [6], which are limit theorem for 
the dynamics of the center of mass of the soliton moving in iV > 2, respectively 
N > 3, under the influence of a homogeneous, time-independent and strongly 
mixing potential satisfying, in addition to assumption (Bl), assumption (B2), 
Section [2J The main ingredient of our analysis is proving (strong) convergence of 
the center of mass dynamics to the auxiliary dynamics of a classical particle in the 
weak-coupling/space-adiabatic limit, LemmattaE] and [7J below. This allows us to 
use the results of [13] . [15] and [TF] for the limiting dynamics of a classical point 
particle in a random potential. For the convenience of the reader, we summarize 
the main conceptual ideas in the analysis of [H], [15] and [17] in the following 
subsection. 
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5.1. Proof of Theorems [5] and [61 We know from TheoremHl that up to times 
of order 0(\logh\/Xh), the soliton center of mass behaves like a point particle 
over spatial and temporal scales of order 0(/i _1 ). This motivates introducing the 
scaling 

a := ha 
v := v 
t := ht. 

It follows from Theorem H] that the rescaled dynamics of the center of mass of 
the soliton, for t G [0, Ce\ log/i|/A), is given by 

(66) c\a = v + 0(h 2 ~ e ) 

(67) = -2XVV(a;uj) + 0(/i 1 ~ e ), 
with initial condition satisfying 

(68) ||a(0) - ha \\ = 0{h 2 ) 

(69) l|v(0)-v o || =OQi). 

We introduce the auxiliary dynamics corresponding to a classical particle in 
the random potential, 

(70) (ha = v 

(71) c\w = -2AVF(a;u;), 
with initial condition 

(72) a(0)=0, v(0)=v . 

This dynamics is "close" to the effective dynamics of the solitary wave in the 
limit h — *■ 0. Moreover, it has been studied extensively as a model of stochastic 
acceleration of classical particles in a random potential, |14j-[17j. 

We note that momentum diffusion occurs at scales of order 0(A~ 2 ), while 
spatial diffusion occurs at scales of order 0(A _2_/3 ), for (3 > 0. This motivates 
comparing the limiting behavior of the true dynamics of the center of mass of 
the soliton to the auxiliary one corresponding to a classical point particle in the 
external potential. We have the following lemma about the convergence of the 
effective dynamics of the center of the soliton to the auxiliary dynamics over 
scales 0(A~ 2 ). 

Lemma 5. Assume (Al)-(A5), (Bl) and (B2), and suppose that there exists 
a > such that A — > as h — ► with \ logh\\ 3 ^ 2+a — >■ oo. Then, for any finite 
T > 0, the stochastic process 

(A 2 a(t/A 2 ),v(t/A 2 )) fe[0 , r] 

converges F-a.s. (strongly) to the stochastic process 

(A 2 a(t/A 2 ),v(t/A 2 )) Ie[0 , T] , 
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as X,h — ► 0. 

Proof. We fix finite T > 0, and choose A and h small enough such that T/X 2 < 
Ce\ log/i|/A. We define 

f Xih (t) :=A 2 ||a(t/A 2 )-a(t/A 2 )||, 

and 

^(t) = ||v(t/A 2 )-v(t/A 2 )||, 
for t <G [0, T]. We have from ([66]) - (E2D that 

(73) \d-Jx,h(t)\ <C(g x , h (t) + h 2 - e ) 
and 

\d- t gx, h \ < C(\\mt/X 2 ) - a(t/A 2 )|| sup || V 2 V\\ L ~ + h^X' 2 ) 

(74) ^^^/a^ + ^A- 2 ), 

where C is a positive constant that is independent of h, e and A. We now use 
Gronwall's lemma and Duhamel formula to obtain bounds on f Xj h and g\,h- 
Introduce the auxiliary C 1 (M) functions fx,h and g x ,h that satisfy 

(75) \dtfx, h (t)\ < C~gxA~t) 

(76) |^ M (t)|<C^/ AA (t), 

with initial conditions fx,h(fy — fx,h(0) and gx,h(0) = gx,h(0). We also introduce 
the rescaled C 1 (1R) functions / and g defined by 

(77) J XA (t) ■= A- 3/2 A,,(A 3 / 2 t) 

(78) g x>h (t):=g x ,h(X 3/2 t). 
They satisfy the differential inequalities 

\dtf Xih (t)\<Cg Xth (t) 
\d,g XA (t)\<Cj Xth (t). 

We define 

h,h '■= f Xl h + 9x,h- 
It follows from the above inequalities that 

\dth,h\ < 2Ch,h- 

By Gronwall's lemma, we have 



sup l x>h (t) < e 2CT l x , h (0), 
te[o,T] 
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which implies 

sup (|7 M (*)i> \9x lh m < e CT (\fx, h (0)\ + 1^(0)1). 

te[o,T] 



Using ( 1781) . we have that 



(79) sup (\f\ h {t)\,\~g x , h m <C'e CT ^ h, 

te[o,T] 

for some constant C that is independent of h, e and A. Now, it follows from 
( T73>fl7~"~) . (E~S) and the Duhamel formula that 

(80) sup (\fx,h(t)l \9x,h(t)\) < CV^'V^A" 2 . 
te[o,T] 

Since | \ogh\\ 3 ^ 2+a — > oo as A, /i — > for some <5 > 0, we have that 

(81) lim e CT ' x3/2 h l ^\- 2 = 0, 

A,/i->0 



which, together with fl80|) imply the claim of the lemma. □ 

This lemma allows us to apply the results of [H] and [15] on momentum dif- 
fusion for weakly random Hamiltonian systems. 

Proof of Theorem^ We know from p~~f] that in dimensions N > 3, and from [15] 
in N = 2, that the auxiliary stochastic process 

(A 2 a(t/A 2 ),v(t/A 2 )) ? > 

converges in law, as A — > 0, to the stochastic process (J * v(s)ds, v(t)), where 
v satisfy (T2TT) . The claim follows now from Lemma [5] and the fact that P-a.s. 
convergence implies convergence in law. □ 

We now specialize to the case N > 3. We introduce the auxiliary Liouville 
equation corresponding to a classical particle moving in the random potential. 
Let A (x, i, k) satisfy the Liouville equation 

(82) d t 4> x = <%% =x ,v=k- V x A + ^v|a= x ,v=k- V k A = k- V X A -2AV X F- V k A . 

with initial condition A (x, 0,k) = o (A 2+/3 x, k), j3 > 0, where 0o is compactly 
supported, twice differentiable in x e M. N , and four times differentiable in k e R , 
such that its support is contained in the shell 

A(M) := {(x, k) G R 2N , 1/M < ||k|| < M}, 

for some M > 1. We now show that the solutions of (1231) and (18*211 are almost 
surely in C\R; C^R™) fl W" 1 ' 00 ^ 2 ^)). 

Lemma 6. Let </> A and A fre solutions of the Liouville equations $23\) and ( ET 
respectively, with the same initial condition <po as described above. Then 

A and A GC 1 (K;C 1 (R M )n'y 1 ' 00 (K M )) P-a.s., 
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for all t > 0. 

Proof. Consider first the Liouville equation corresponding to the auxiliary dy- 
namics, (|82[) . The analysis for f[2"3"j) is similar. Introduce the variables X := 
(x,k) G R 27V , and let A(X) = (k, -2AV x F(x))| X: =(x,k), $(X,t) := A (x,i,k). 
Note that 

A(X) = A l (X)+A n (X), 

where A { is linear in X and A n is nonlinear in X, with A n G C 1 (R 2JV )nW' 1,00 (]R 2JV ) 
as a vector- valued function, since V G W 2,QC H C 2 . Now, Eq. fl82l can be written 
as the Hamilton- Jacobi equation 

^$(X,t) = A(X)-V x <f(X,t), 

with initial condition $ £ C 2 (M 2Ar ). We use the method of characteristics to 
solve the above Hamilton- Jacobi equation; see for example [31] for a discussion 
about Hamilton- Jacobi equations. Consider the mapping 

Xt,*o(X):=X-A(X)(t-to) eK 2 " 

We want to show that P-a.s. x is a diffeomorphism of class C l , i.e., for w 6 fi, 

Y is biiective and with a C 1 inverse. For u; G f2 and < t — to < 7; hnni — • 

the mapping \ is invertible and differentiable, since the linear part A\ is clearly 
invertible and differentiable, while the nonlinear part A n is a differentiable per- 
turbation for short times. 

We claim that for any t — to > 0, the mapping x is a P-a.s. diffeomorphism 
of class C 1 . Let r > be the maximal time such that Xr+t ,t is P-a.s. invertible 

with C l inverse, and suppose that r < oo. Let St := -. hnni — • It follows 

from the definition of x that 

.St St 
X T +%,t = Xr-f ,t + Ay + 

Since r — ^ < r, x r -^i i s a ^ diffeomorphism. Furthermore, the linear part 
is differentiable and invertible. Again, A n ^ is a small differentiable perturba- 
tion, and hence X T +^ t is a diffeomorphism. However, this contradicts the 
definition of r, and hence r = oo. 

Using the definition of x> h is straightforward to verify that the solution $ of 
the Hamilton- Jacobi equation satisfies 

$(x t ,o(X),t) = $ (X). 

Since x has a differentiable inverse P-a.s., 

$(x,t) = $ (x,; 1 (x)). 

Furthermore, $ G C 2 (M 27V ), implies that $ G C X (R; (^(M 2 ^) n W 1>0O (R 2JV )) P- 
almost surely. 
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Using the same argument, one can show that the solution of fl23|) with initial 
condition 4> is almost surely in C X (E; C^IR 2 ^) n W 1,O0 (R 2N )) by studying the 
solution of the corresponding Hamilton- Jacobi equation. □ 

We also have the following lemma on the convergence of the (rescaled) solutions 
of ( 123]) and ( 1821) in the weak-coupling/space-adiabatic limit. 

Lemma 7. Let <p x and <p x be solutions of the Liouville equations and (EIP, 
respectively, with the same initial condition (f> as described above. Suppose that 
there exists a > such that A —>■ as h — > with \ logh\X 1+a — > oo. Then, for 
any fixed T > and any j3 G (0, a/2), we have that 
(83) 

lim sup [0 A (x/A 2+/3 ,t/A 2+2/3 ,k)-0 A (x/A 2+/3 ,t/A 2+2/3 ,k)] = O P-a.s.. 

A >^0 (i,x,k)£[0,T]xA' 

Proof. For fixed < T < oo, let A, /i be small enough such that TA -2 ~ 2/3 < 
Ce\ \ogh\/X. This is possible since | \ogh\X 1+a — ► oo as A, /i — > 0, and /3 < a/2. 
Recall that A and cf) x satisfy the same initial conditions. Let 




It follows from (|23|) . (|82|) and Lemma [6] that, 

<9^ A = k ■ Vj X - 2AV x F(x) ■ V k A + 0(^11^11^^)), 

for t G [0,T], with initial condition A (O) = 0. Therefore, for t G [0, T] and all 

uj E Q, 

|0 A (x/A 2+/3 ,t/A 2+2/3 ,k)| = 0(TA- 4 - 3/ V- e ). 
Now, for < (3 < a/2, 

lim T\- A ~ w h l - e = 0, 

which implies the claim of the lemma. □ 

We now use Lemma H and the results of p2] to prove Theorem [6j 

Proof of Theorem It is shown in jTTJ that momentum diffusion of a classical 
particle in a random and strongly mixing potential converges to a spatial Brow- 
nian motion over longer time scales (see Subsection 15.21) . As a cosequence, the 
auxiliary stochastic process given in (|7Dl)- (172l) converges in the weak-coupling 
limit to a spatial Brownian motion: For A<1, there exist (3 G (0, a/2) and some 
constant C independent of h and A, such that, for all < /3 < (3, fixed <t < T 
and all compact sets K C A(M), we have 

(84) sup \E[4> x (x/\ 2+ P, t/A 2+2/3 , k)] - it(x, t, k) | < CTX^" 13 , 

(t,x.,\i)£[to,T]xK 
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where A satisfied (1821 and u satisfies (124|) . We also have from Lemma [7] that 

(85) lim sup E[0 A (x/A 2+/3 , t/A 2+2/3 , k) - A (x/A 2+/3 , t/A 2+2/3 , k)] = 0. 
M-*o (t )X ,k)e[t ,r]xK 

It follows (ED and (EH} that 

lim sup |E[</> A (x/A 2+/3 , t/A 2+2/3 , k)] - u(x, t, k) | = 0. 

A ' fe ^ (t,x,k)6[t ,T]xl<: 

□ 

5.2. Remarks on diffusion for weakly random Hamiltonian flows. In this 
subsection, we recall the main conceptual ideas of references [H], [15] and [17j . 

Let (a A (t), v A (t)) be the rescaled positions and velocities of a classical particle 
in a random Hamiltonian flow defined by 

a A (t) := A 2 a(t/A 2 ), 

v A (t) :=v(t/A 2 ). 

They satisfy the differential equations 

<%a A = v\ 

<%v A = -^VF(a A /A 2 ), 

with initial condition 

a A (0) = 
v A (0) = v . 

The main difficulty in obtaining the limit of the above stochastic process is that 
a A (t + St), Si <C 1, may be correlated to a A (t), and hence the process may be 
non- Mar kovian . 

We start with discussing the case iV > 3, which was studied in [TJ]. The 
authors of [2] introduce an auxiliary modified dynamics and a stopping time 
r\ such that, up to times T\, the momenta of the modified process are locally 
aligned and the modified trajectory is a straight line during times of intersection. 
Furthermore, for times larger than t\, the modified process is a diffusion process. 
Because of the spatial mixing property of the random potential, the limit of the 
modified process is Markovian. The proof is completed by showing two more 
elements: First, the process (a A (t), v A (t)) 0< j <TA converges weakly to the modified 
process as A — > 0. Second, the stopping time r x — > oo as A — > 0. 

The case of N = 2 is a little bit more difficult, due to the fact that the 
limiting process self-intersects. This difficulty is overcome in [TJ] by modifying 
the stopping time condition for the modified process: up to stopping time t\, 
only transversal intersections of the modified trajectory are allowed. 

We also recall the main ingredients in [T7] for proving convergence of the tra- 
jectory of a classical particle in a strongly mixing random potential in dimensions 
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iV > 3 to a spatial Brownian motion in the weak-coupling limit. The proof is 
based on extending the analysis in [T3] to obtain explicit estimates on the con- 
vergence to momentum diffusion, and then using asymptotic expansion to obtain 
spatial diffusion from momentum diffusion on longer time scales. 
Consider e Cl' h2 (R N x [0, +oo) x M w \{0}) satisfying 

JV 

(86) d t J = d h Dij(k)d k .0 + k • V x 0, 

with initial condition 0(x, 0, k) = 0o( x ; k). Using an extension of the analysis in 
[H] and [30], it is shown in [17] that there exists (3 > such that, for all compact 
sets K e A(M), 

(87) sup |E[0 A (x/A 2 ,t/A 2 ,k)] -0(x,t,k)| < CT(1 + \\Ml^) X ^ 

(t,x,k)e[0,T]xiC 

where C is independent of A and T. Note that (15B1) is well-posed in C b ' ' 2 , see for 
example [32] . Furthermore, the nonvanishing property of the correlation function 
implies that the diffusion matrix has rank N — 1, and hence (1861) corresponds to 
diffusion on a sphere of constant momentum, see [H] and [TTJ. By applying stan- 
dard asymptotic expansion techniques, one can show that the long-time limit of 
the solution of flHBI) is spatial diffusion: For every < t < T < oo, 0(x/7, t/7 2 , k) 
converges in C([t ,T\; L°°(R N xR N )) as 7 ^0 to u(x,t,k), where u satisfies (T2Ij). 
with 

(88) sup | M (x,t,k)-0(x/ 7 ,t/ 7 2 ,k)| <C( 7 T + ^7)110011^, 
te[to,r| ) (x,k)eR-tfxR JV ' 

where C is a constant that is independent of A and T. Now, (IHTj) and (|88|) yield 

(HU). 

6. Appendix A: Well-posedness of a generalized nonautonomous 

NONLINEAR SCHRODINGER EQUATION 

We now discuss the local and global well-posedness of a generalized nonau- 
tonomous NLS equation with random (time-dependent) nonlinearities and po- 
tential. Global well-posedness and possible occurence of blow up are addressed 
in [33] and [31] for the NLS equation with power nonlinearities and additive or 
multiplicative random potential in the form of white noise. 

Consider the problem corresponding to a generalized nonlinear Schrodinger 
equation 

(89) idtip = -Aip + g(t,if);w),ip(t = 0) = <l>, 

where u G Q and g contains both the potential and the nonlinearity. Here, g 
also depends on x e Mr, but we drop the explicit dependence when there is no 
danger of confusion. 
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We say that (q, r) is an admissible pair if 

2N 

rG [ 2 'iV^)> (re[2,oo],lV = l) 

(90) \ = 

We make the following assumptions on g. 
(CI) The nonlinearity g = g 1 + ■ ■ ■ + with 

gj G C(R, C(H\ H' 1 )) P - a.s., j = 1, • • • , k. 

(C2) There exist admissible pairs ((ft, r,-), j = l,---k, such that, for every 
T, M > 0, there exist a constant C(M) independent of T, and (3 indepen- 
dent of T and M, such that 

||<&(*,«)-<&(M)II £ ^ W) <C(M)(l + T^\\u-v\\ L r HRN) , P-a.s., 

for all u,v E H 1 with ||tt||^i + II^Hh 1 < ^ an d |£| < T", where r' is the 
conjugate of r, i.e., 1/r + 1/r' — 1. Furthermore, 

ll&MII^ <C(M)(l + ^)(l + |H| w i„,), P-a.s., 

for all u G i? 1 (R JV ) fl W 1 ' r (R JV ) such that \\u\\ H i < M and \t\ < T. 

(C3) 

Img^t, w)w = 0, j = 1, • • • , fc, P — a.s. 

almost everywhere on R-^, for all t G R and u G if 1 . 
(C4) There exists a functional G, G C(R, C^if^R)) with G$ = where the 
prime stands for the Frechet derivative. We let G — G± + • • • G^. For 
u G H\ 

(91) \d t G(t,u;uj)\ <C(\\u\\ L 2)l(t)¥-a.s., 

where C depends only on ||m||l2 and the real function I G L°°(M) such 
that l(t) < 1 for almost all tGl. 
(C5) For all M > 0, there exists C(M) > and e G (0, 1), both independent of 
t G R, such that 

(92) |G(*,u;w)| < -^\\u\\ 2 H i+C(M), P-a.s., 

uniformly in t G R, Vu G if 1 , such that \\u\\ L 2 < M. 

In what follows, we let Ocfi, with fx(Q) = 1, denote the set over which (Cl)- 
(C5) hold. We have the following result about local well-posedness almost surely 
in H 1 . 

Proposition 2. Suppose g satisfies assumptions (C1)-(C3). Then the following 
holds. 
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(i) For every <ft G H 1 (M. N ), there exists F-a.s. strong H 1 -solution u of / fgpj) . 
which is defined on a maximal time interval (— T*,T*), such that there 
exists a blow-up alternative, i.e., if T* < oo, ess — sup tJgn ||M(t)||//i — > oo 
as t y T*, and if T* < oo, ess — sup wgn ||w(t)||#i — > oo as t \ — T*. 
Moreover, 

/or a// admissible pairs (a, 6). 

(ii) The charge is F-a.s. conserved, 

II u WIIl2 = \\<f>\\v», a; G Q, 
for allte (-%,T*). 

(iii) F-a.s., u depends continuously on <fi : If 4> n ™— >° 4> ^ n H 1 , and if u n is the 
maximal solution of / fffPj) corresponding to the initial condition <fr n , then 
u n n —>° u in C([—S,T],L P (M N )) for every compact interval [—S,T] C 
(-T*,T*) andpe[2,^- 2 ) (p G [2, oo), N = 1). 

The proof of the above proposition is a direct extension of the deterministic 
case with u G f2; see the Appendix in [7] for a discussion of the latter. It is 
based on Kato's method, which relies on Strichartz estimates and a fixed point 
argument, [25l 126]. 

Proving global well-posedness for data which are not necessarily small is a 
little bit more delicate, since energy is not conserved. We now define the energy 
functional 

(93) E(t, u; u):=^J \Wu\ 2 dx + G{t, u; lu), 

for u G if 1 (M Ar ) and uj G fl. Since the nonlinearity and the potential depend on 
time, the energy is not conserved. We have the following proposition. 

Proposition 3. Suppose that (C1)-(C4) are satisfied, and let u denote the solu- 
tion of / fggj) given by Proposition Then 

(94) \E(t,u(t);u;)\ < \E(0,fru)\ +TC(\\<f>\\ L2 ), P-a.s., 

for allt G [— T, T], where [— T, T] is a compact subset of (— T*, T*), and C( 1 1 1 1 z>a ) 
appears in (C4). 

Proof. Let to G Q. Since (C1)-(C3) are satisfied, the results of Proposition [2] hold. 
We choose a finite T > such that T < mh^T*, T*). We have from Proposition 
[2] that 

M GLL((-T„T*),iy 1 ' b (M 7V )), W eH 
for all admissible pairs (a, 6) . In particular, 

u G L qj ( (— T, T) , W 1,r J (R^) ) , j = 1, • ■ ■ , fc, wGfi 

where the admissible pairs (qj,Tj) appear in (C2). 
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Since V commutes with the unitary propagator U the L 2 norm is invariant 
under unitary transformations, we have using the Duhamel formula that 

\\Vu(t)\\ 2 L2 = \\WU(0,t)u(t)\\ 2 L2 



t 

2 



dsU(0,s)Vg(s,u(s);uj)\\ 2 L 2 

Jo 

\V(j)\\ 2 L2 - 2Im(V0, / ds U(0,s)Vg(s,u(s);u)) + \\ [ ds U(0 } s)Vg(s,u(s);uj)\\ 2 L2 
Jo Jo 

|V0||£ 2 + 2Im/ ds (V^(a,«(a);w),t/(s,O)V0> + || / ds U(0, s)Vg(s, it(s); w) ||£ 2 
Jo Jo 

|V0||i 2 + 2Vlm / (V ft (s,u(s); W ),VM(s)) 
i=i " /o 
fc /■* 

|V0||^ 2 - 2 V"lm / {gj(s,u(s);uj), Au(s)), cu eU 

where the scalar product is well-defined using (C2) and duality on 

(L\(-T, T), H l )+L^{{-T, T), H 1 ^)) x (L°°((-T, T), ff 1 )nL*((-T, T), H 1,rj ))), 

j — 1, ■ ■ ■ ,k, see for example [29]. Here, 

H S ' P (R N ) := {u G : JF^l + \k\ 2 )^u G ^(R^)}, s G R, 1 < p < oo, 

where T stands for the Fourier transform; and the space H s,p is equipped with 
the norm 

\\u\\ H s, P = + \\k\\ 2 )^u\\ LP , u G H S,P (R N ). 

Now, 

lm(g(t,u(t);uj),Au(t)) = limlm((l - eA) _1 ^(t, u(t); u), (1 - eA) _1 Aw(t)) 

= lim Im((l - eA)-^(t, u(t)), (1 - eA^H^t) + «(*)))) 

= Re{g(t,u(t);u),d t u(t)) 

= ^-G(t,u(t);u) - (d t G)(t,u(t);uj), to G H 
at 

for almost all i G (— T, T), where G appears in Assumption (C4). Therefore, 
(95) 

||Vu(t)||| 2 = ||V0||i 2 -2G(t, M (t);^) + 2G(O,0;u;) + 2 / ds (d s G)(s,u(s);u), 

Jo 

for u G ft. Together with ( |93j) . Assumption (C4) and conservation of charge, this 
implies 

\E(t, u(t),uj)\ < \E(0, <f>;u)\ + TCiUWv), 
for t G [-T, T] and all oo G ft. □ 
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We have the following theorem about P-a.s. global well-posedness in H l . 

Theorem 7. Suppose (C1)-(C5) hold. Then the solution u of [8^1 with initial 
condition (f> G H 1 (M N ) is P-a.s. global in H 1 , i.e., foruj G Q, 

T* = T* = oo, 

where T*,T* appear in Proposition^ 

Proof. Let uj G Q. Propositions [2] and [3] follow from (C1)-(C4). It suffices to show 
that ||w(£)||#i, t G [0,T*) is finite if T* < oo, which, together with the blow-up 
alternative, implies a contradiction. Suppose T* < oo. Assumption (C5), charge 
conservation and flM|) imply that 

^IK«)II 2 ^ = ^(IK«)II! 2 + I|v^)|| 2 i2 ) 

< l\Hs)\\ 2 L 2 + \E(s,u(s);u))\ + \G(s,u(s);u))\ 

< \u\\b + \mM\ + t*c(U\\l*) + ^Hs)\\ 2 Hl + c(u\\ L2 ), 

for all s G [0,T*), ufSl, and hence 

sup sup Hh 1 < oo, 

a, e ns6[0,T*) 

for finite T*, which contradicts the blow-up alternative. The case of T* is proven 
similarly. 

□ 

7. Appendix B: Rate of change of field energy and momenta 

Proof of ( E§P ; Section^ Differentiating (if), X Vhif)) with respect to t and using 
dl]), Assumption (Al), and the fact that if) G i/ 1 (lR iV ), we have, for uj G Q, 

d t (if), XV h if)) = (if), Xd t V h iP) + (d t if), \V h ifj) + (if), \V h d t if)) 

= (if), Xd t V h iP) + (-Aif) + \V h if) - f{if)),i\V h if)) 
+ (iif),\V h (-Aif; + \V h if)-f(iP))) 
= (if;,Xd t V h if)) + 2(iXVV h if;, Vif)), 

where we have used integration by parts in the last step. □ 

Proof of ( f57y ; Section [7} The proof of floT|) follows directly from flUB^) in the proof 
of Proposition [3j Appendix A, with the identification 

g(t,u) = XV h (t)u- f(u). 

□ 



:;<s 



Proof of (EE) , Section [7} We use a regularization scheme similar to the one used 
in Proposition [3j Appendix A. Let I e := (1 — eA) _1 . For all u G Q, we have that 

d t (iil), Vip) = d t lim(/ e z?/>, I e Vip) 

e\0 

= \im{(I e id t ip, J e VV) + (Icfy, IeVd t ^}} 

= lim{(/ e (-AV> + AW - f(ip)),I e ViI>) - (J e V,/ e V(-A^ + AW - 
= A(V,V r h VV)-A(V,V(W)) 

= -%v^). 

□ 
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